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ABSTRACT 

We present analytic solutions of Maxwell equations in the internal and external background 
spacetime of a slowly rotating magnetized neutron star. The star is considered isolated and in 
vacuum, with a dipolar magnetic field not aligned with the axis of rotation. With respect to a 
flat spacetime solution, general relativity introduces corrections related both to the monopolar 
and the dipolar parts of the gravitational field. In particular, we show that in the case of infinite 
electrical conductivity general relativistic corrections due to the dragging of reference frames 
are present, but only in the expression for the electric field. In the case of finite electrical 
conductivity, however, corrections due both to the spacetime curvature and to the dragging of 
reference frames are shown to be present in the induction equation. These corrections could 
be relevant for the evolution of the magnetic fields of pulsars and magnetars. The solutions 
found, while obtained through some simplifying assumption, reflect a rather general physical 
configuration and could therefore be used in a variety of astrophysical situations. 

Key words: relativity - (magnetohydrodynamics) MHD - stars: neutron - rotation - mag- 
netic fields 



1 INTRODUCTION 

The investigation of the influence of strongly curved spacetimes on the properties of electromagnetic fields has an interest of its own which is 
further increased when these effects could be related to a rich observable phenomenology. This coupling between general relativistic effects 
and electromagnetic fields is expected to be particularly important in the vicinity of neutron stars which are among the most relativistic astro- 
physical objects and are characterized by very intense magnetic fields (Lamb 1991, Glendenning 1996). A number of different observations 
indicate that in young neutron stars the surface magnetic field strengths are of the order of 10 11 - 10 13 G. In some exceptional cases, as 
those of magnetars, magnetic field strengths > 5 x 10 14 G are considered responsible for the phenomenology observed in soft gamma-ray 
repeaters (Duncan & Thompson 1992, Thompson & Duncan 1995). Older neutron stars, observed as recycled pulsars and low mass X-ray 
binaries, show instead surface magnetic fields that are much weaker < 10 10 G suggesting that these are subject to a decay, even if it is still 
difficult to establish whether the decay is due to accretion (Geppert & Urpin, 1994; Konar & Bhattacharya, 1997) or to other processes. 

In the case of isolated neutron stars, the possibility of magnetic field decay as a result of accretion does not arise, but there are still a 
number of different ways in which the energy stored in the magnetic can be lost. This can happen either through the emission of electromag- 
netic (dipole) radiation, through Ohmic decay, through ambipolar diffusion, or through more complicated effects such as "Hall cascades" 
(see Goldreich and Reisenneger 1992 for a review). The investigation of these scenarios requires combined efforts. On one hand, there is 
the search for a more precise description of the microphysics of the processes involved, some of which are still not well quantified. On the 
other hand, attention is paid to a more realistic description of the gravitational effects on the properties of the electromagnetic fields in highly 
curved spacetimes and this is also the motivation of this work. 

The investigation of the general relativistic corrections to the solution of Maxwell equations in the spacetime of a relativistic star has 
a long history. The initial works of Ginzburg & Ozernoy (1964), Anderson & Cohen (1970) and of Petterson (1974) on the stationary elec- 
tromagnetic fields in a Schwarzschild spacetime have revealed that the spacetime curvature produces magnetic fields which are generally 
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stronger than their Newtonian counterparts (see also Wasserman & Shapiro 1983 for a subsequent derivation). Sengupta (1995) has reconsid- 
ered this problem and also looked for a general relativistic expression for the electric field in the Schwarzschild background of a neutron star. 
As we will discuss in Section 3.2.2 the method used in his derivation is not entirely correct and the results obtained for the electric field are 
not solutions of Maxwell equations. More recently, Sengupta has also considered the problem of the Ohmic decay rate in a Schwarzschild 
spacetime (Sengupta, 1997). His approach is strictly valid only for the region of spacetime external to the star as it does not provide a correct 
general relativistic description of the electromagnetic fields internal to the star. Within these approximations, however, Sengupta (1997) has 
pointed out that the effects of intense gravitational field seem to decrease the overall decay rate by a couple of orders of magnitude. The 
same problem has also been considered in more detail by Geppert, Page and Zannias (2000). Their analysis was aimed at a mathematically 
consistent solution of Maxwell equations also in the spacetime region internal to the star and makes therefore use of a generic metric for a 
non-rotating relativistic star. Their results, while confirming a decrease in the typical decay time for the magnetic field, also show that the 
decay time is smaller but comparable with the one found in flat spacetime. 

The general relativistic effects induced by the rotation of the star were first investigated by Muslimov & Tsygan (1992) in the slow 
rotation approximation. A similar approach was also used by Muslimov and Harding (1997) for the electromagnetic fields external to a 
rotating magnetized star. Their analysis refers to a charge filled magnetosphere and represents the relativistic extension of the Goldreich- 
Julian model. Using a different derivation, Prasanna and Gupta (1997) have also investigated the properties of the electromagnetic fields in 
the magnetosphere of a relativistic rotating neutron star, with special attention being paid to the dynamics of charged test particles. 

We here extend and unify all of the above investigations by considering the solution of Maxwell equations in the internal and external 
background spacetime of a slowly rotating magnetized relativistic star. The star is considered isolated and in vacuum, with a dipolar magnetic 
field which is not assumed aligned with the axis of rotation. The purpose of this paper is threefold. Firstly, we want to extend previous results 
to the most general case of a misaligned rotator, providing for this case also the form of the electric field. Secondly, we want to discuss the 
possible role played by frame dragging effects in the Ohmic decay for an isolated neutron star and estimate its importance. Thirdly, we wish 
to clarify a few important aspects of the solution of Maxwell equations in the gravitational field of a relativistic star that, when overlooked, 
have led to incorrect solutions (Sengupta 1995, Prasanna and Gupta 1997). Finally, by providing a rather general solution to the problem 
(although truncated at the lowest order in the expansion of the angular dependence) we offer a compact reference from which all of the 
previous results can be easily found in the appropriate limits and which could have practical astrophysical applications. 

The paper is organized as follows: in Section 2 we write the general relativistic Maxwell equations in the metric of a slowly rotating star 
and the form they assume when the electromagnetic fields are those measured in the orthonormal frame of zero angular momentum observers. 
In Section 3 we find the stationary solutions (i.e. solutions in which the infinite conductivity of the medium prevents a variation in time of 
the star's magnetic moment) to Maxwell equations outside and inside the misaligned rotating star. For this we consider first the problem in 
Newtonian gravity and we then extend the results to general relativity within the slow rotation approximation. Section 4 is devoted to the 
equivalent problem, but in the case in which the magnetic field is not supposed stationary. There, we derive the basic induction equations for 
the evolution of the inner stellar magnetic field of a misaligned rotating star. Section 5 contains our conclusions and the prospects of future 
developments. 

A number of appendices provide further details about some of the calculations carried out in the main part of the paper. In particular, 
Appendix A summarizes the components of the electromagnetic tensor in a coordinate basis and in a locally orthonormal tetrad, while 
Appendix B shows the derivation of the radial eigenfunctions for the electromagnetic fields in terms of Legendre's equation. Appendix C 
shows the explicit expressions for the surface charges and currents and, finally, Appendix D contains an alternative and equivalent derivation 
of the equations for the time evolution of magnetic field in terms a vector potential. Throughout, we use a space-like signature (—,+,+,+) 
and a system of units in which G = 1 = c (However, for those expressions with an astrophysical application we have written the speed of 
light explicitely.). Greek indices are taken to run from to 3 and Latin indices from 1 to 3; covariant derivatives are denoted with a semi-colon 
and partial derivatives with a comma. 



2 MAXWELL EQUATIONS IN A SLOWLY ROTATING SPACETIME 

The difficulties of an analytic solution of the Einstein-Maxwell equations in the proximity of a rotating relativistic star inevitably force us to 
the use of some approximations. The first approximation comes from neglecting the influence of the electromagnetic field on the metric and 
by solving Maxwell equations on a given, fixed background*. The second approximation is in the specific form of the background metric 
which we choose to be that of a stationary, axially symmetric system truncated at the first order in the angular velocity Q. In a coordinate 
system (ct, r, 6, cf>), the "slow rotation metric" for a rotating relativistic star is (see, for example, Hartle 1967, Hartle & Thorne 1968, Landau 
&Lifshitz 1971) 

ds 2 = -e 2 * (r) di 2 + e 2A(r) dr 2 - 2u(r)r 2 sin 2 6dtd<f> + r 2 d8 2 + r 2 sin 2 6d<j} 2 , (1) 



* This is indeed a very good approximation since even for very highly magnetic neutron stars the electromagnetic energy density is much smaller than the 
gravitational one. 
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where uj(r) can be interpreted as the angular velocity of a free falling (inertial) frame and is also known as the Lense-Thirring angular 
velocity. The radial dependence of u> in the region of spacetime internal to the star has to be found as the solution of the differential equation 

1 d ( 4 -duj\ ^dj 
r 3 dr V ^ dr J dr- L 



Zs-Jr]-)+4^ = 0, (2) 



where we have defined 

I = e-(*+ A \ (3) 

and where 

Q = £l-w, (4) 

is the angular velocity of the fluid as measured from the local free falling (inertial) frame. In the vacuum region of spacetime external to the 
star, on the other hand, tu{r) is given by the simple algebraic expression 

goj> _ 2J 
dt g^ r 3 

where J — I(M, R)Q, is the total angular momentum of metric source as measured from infinity and I(M, R) its momentum of inertia (see 
Miller 1977 for a discussion of / and its numerical calculation). Outside the star, the metric (1) is completely known and explicit expressions 
for the other metric functions are given by 

e2Hr) ^ / _ 2M\ = e _ 2A(r) ^ r>R ^ (6) 



w(r)s ^ = _*» " (5) 



r 

where M and R are the mass and radius of the star as measured from infinity. 

An important aspect, often overlooked in the literature, should now be underlined. The metric (1) is the simplest metric that provides all 
of the most important general relativistic corrections to the solution of the Maxwell equations in the gravitational field of a rotating relativistic 
star. The use of a Schwarzschild metric in place of (1) (Sengupta 1995, 1997) is potentially very dangerous. Firstly, and as pointed out by 
Geppert et al. (2000), a Schwarzschild metric allows for a proper treatment of the electromagnetic fields only in the spacetime region external 
to the star and leaves unsolved the problem of a matching of the external electromagnetic fields with the internal ones. Secondly, and despite 
different claims (Sengupta 1997), a Schwarzschild metric is intrinsically inadequate to describe physical systems such as pulsars in which 
the coupling of electromagnetic fields and rotation is a key feature. Note, on the other hand, that using the slow-rotation approximation gives 
rather accurate results for all pulsar periods so far observed. The metric (1) has coefficients each of which is the lowest-order term of a series 
expansion in ascending powers of 57. Comparing the magnitude of the neglected higher order terms with that of the one retained in each case, 
gives ratios of the order R 3 Q 2 /GM which is smaller than 10% even for the fastest-known millisecond pulsar PSR 1937+214. 

The general form of the first pair of general relativistic Maxwell equations is given by 

3!F Mi7] = 2 (F Q(3 , 7 + F 7Qi/3 + F fj -. ha ) = . (7) 

where F a/ 3 is the electromagnetic field tensor expressing the strict connection between the electric and magnetic four- vector fields E a , B a . 
For an observer with four- velocity 11°, the covariant components of the electromagnetic tensor are given by (Lichnerowicz 1967; Ellis 1973) 

F afj = 2u [a E g] + rj al 3~ / 6U~ f B S . (8) 

where TJ Q(3 ] = ^(T Q/ 3 — Tpa) and r/c^s is the pseudo-tensorial expression for the Levi-Civita symbol e Q( 3 7 ,5 (Stephani 1990) 

a8^/5 1 i /n\ 



'-9 

with g = dct\g a fj\ = — e 2 (*+ A ) r 4 s i n 2 Q for the metric (1). A useful class of observers is represented by the "zero angular momentum 
observers" or ZAMOs (Bardeen, Press & Teukolsky 1972). These are observers that are locally stationary (i.e. at fixed values of r and 6) 
but who are "dragged" into rotation with respect to a reference frame fixed with respect to distant observers. At first order in Q, they have 
four- velocity components given by 



( W «) ZAMO =e-*M(l,0,0 )W ) ; KOzamo =e* W (- 1,0,0, o) 



(10) 



In the coordinate system (ct, r, 6, <f>) and with the definition (8) referred to the observers (10), the first pair of Maxwell equations (7) 
take then the form (see Appendix A for the explicit expressions of the electromagnetic tensor) 

(eVsiri6»/r) =0 , (11) 

(eV sin €>)^ = e* - E^g) - (c^eV sin o) B r A , (12) 

(eV sin 0)^= fa e*) - e*E ri4> - (c^eV sin o) B% , (13) 
(e V sin 0) ^ = - (Eb e*) + e*E T ,e + sin B (o;e A r 2 B r ) + V (sin 9B e ) g . (14) 
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The general form of the second pair of Maxwell equations is given by 

F al3 = 4nJ a (15) 
where the four-current J a is a sum of convection and conduction currents 

J a =p e w a +j a , j a w a =0, (16) 

with w being the conductor four-velocity and p e the proper charge density. If the conduction current j a is carried by the electrons t with 
electrical conductivity a, Ohm's law can then be written as 

j a = aFafjw* 3 , (17) 
while a more general expression can be found in Ahmedov (1999). We can now rewrite the second pair of Maxwell equations as 

(eV sin 60S 1 ) = 47re 4>+ V sin 6 J° , (18) 

e *(^, fi -ft^)-( W eVsiii9)£^ = (eVsinfl) ^ + 47re*+V sin0J rr , (19) 

e*B r ,* - (e*B*) - (weV sin d) B% = (eV sinfl) ^ + 47re* + V sin^j" , (20) 

(e*B e ) - e *B r , e + sin^^ e V£ r ) + ueV (sin ^ = (e V sin 6>) + 4jre*+ V sin 6> J* . (21) 

Maxwell equations assume a familiar flat-spacetime form when projected onto a locally orthonormal tetrad. In principle such tetrad is 
arbitrary, but in the case of a relativistic rotating metric source a "natural" choice is offered by the tetrad carried by the ZAMOs. Using (10) 
we find that the components of the tetrad {ep} — (eg, e? , eg, e^) carried by a ZAMO observer are 

e? =e~*( 1,0,0, J) , (22) 



e" = e~ A ( 0. 1.0.0 ) . (23) 
e? = i(0,0,l,0) , (24) 



r 



0,0,0,1 . (25) 



The 1-forms {cj m } = (u> , u> r , u> , u> ), corresponding to this tetrad have instead components 

w° = e* ( 1,0,0, oV (26) 



u£ =e A (0, 1,0,0 ) , (27) 



(28) 

ujI = r sin 6 ( - w, 0, 0, 1 ) . (29) 



We can now rewrite Maxwell equations (1 1)— (14) and (18)— (21) in the ZAMO reference frame by contracting (7) and (15) with (22)-(25) 
and (26)-(29). After some lengthy but straightforward algebra, we obtain Maxwell equations in the more useful form 

sin 9 (r 2 B f ^ + e A r (sin 8B^ + e A rB^ = , (30) 



(r sin 8) = e 



- (sin 603^) 



(wr sin 0)5^, (31) 
(e A r sin fl) ^ = -e* +A S^ + sin (re*£*) - (toe\ sin 6>) S% , (32) 
(eV)^ = -(re*S*) + e* +A £% + sinfl (cr 2 B f ) + u,e A r (sin ^ (33) 

t This is a reasonable assumption if the neutron star has a temperature such that the atomic nuclei are frozen into a lattice and the electrons form a completely 
relativistic, and degenerate gas. 
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and 



sin 9 {fET) + e A r (sin 9E^ ^ + e A rE^ 



{lot sinO) E\ 



e* (sinOB^ g -B f 
e* +A B% - sin 9 (r e*B*) - (ue A r sin (?) E% 
(e*r.B*) - e*+ A B% + sin 9 (cor 2 E^ + ioe A r (sin 9E^ ^ 



4ne A r 2 sin 6 J 1 



dE r 

(r sin 9) — h 4-7re V sin 9 J r 

at 



= (e A r sin fl) + 4^e*+ A r sin 6 J« , 
= (e A r) 



+ 47re* +A r J"* + 47re A cjr 2 sin 0J f 



(34) 
(35) 
(36) 
(37) 



Equations (35)-(37) can now be rewritten in a more convenient form. Taking our conductor to be the star with four-velocity components 



* w ( 1,0,0, n 



*(r) 



1,0,0, 



p2*(r) 



we can use Ohm's law (17) to derive the following explicit components of J a in the ZAMO frame 

j ur sin 9 i 

J = pe + a — -=- — E y , 



g E 



tor sm f § 
r — B 



r = *iE e + Qr T e B < 



(38) 

(39) 
(40) 
(41) 
(42) 



Next, we discuss a few assumptions that are going to be used hereafter. Firstly, we assume there is no matter outside the star so that the 
conductivity a — for r > R and that a ^ only in a shell with R IN < r < R (e.g. the neutron star crust). Secondly, we consider a to be 
uniform within this shell (Note that this might be incorrect in the outermost layers of the neutron star but is a rather good approximation on 
the crust as a whole.). Thirdly, we ignore the contributions coming from displacement currents. The latter could, in principle, be relevant in 
the evolution of the electromagnetic fields, but their effects are negligible on timescales that are long as compared with the electromagnetic 
waves crossing time. In view of this, we will neglect in (35)-(37) all terms involving time derivatives of the electric field and use Ohm's law 
to rewrite equations (35) and (36) as 



r sin 9E r = 



1 

47TCT 



e r sin 9E = 



47TCT 



(sin9B^ g -B S , 
B a sin t 



■O(0) 



+ 0(fi) 



Substituting (43) and (44) in the left hand sides of equations (35)-(37) eliminates the dependence from the electric field and yields 



$+A r>r 



B a. sin ^ 



(sin9B^ ^ 



B 



4na 



B" 



4lTO 



e B sin 



(sin 

(e*rB*) 



(eW)^e- A ^ + i-{^[(sin^)^<]} / 



4-ko 



e i+A B r rl , — sin i 



(e*rB*) 



47rrsin#J r , 
47re* +A rsine/ , 



47re <E>+A r + 47ve A cur 2 sin9p e 



(43) 
(44) 

(45) 
(46) 

(47) 



3 STATIONARY SOLUTIONS TO MAXWELL EQUATIONS 

In this Section we will look for stationary solutions of the Maxwell equation, i.e. for solutions in which we assume that the magnetic 
moment of the magnetic star does not vary in time as a result of the infinite conductivity of the medium. Note that this does not mean the 
electromagnetic fields are independent of time. As a result of the misalignment between the magnetic dipole \i and the angular velocity 
vector 1~2, in fact, both the magnetic and the electric fields will posses a periodic time dependence produced by the precession of pi around i~i 
(see Fig. 1). 
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Figure 1. Schematic representation of a misaligned rotator. Here (eg, e f , eg, e^) is a local orthonormal frame, fj, is the magnetic dipole moment of the star, 
X is the inclination angle relative to the rotation axis, and A the instantaneous azimuthal position. 



3.1 Rotating Magnetized Conductor in a Minkowski Spacetime 

Before looking at the problem of a magnetized rotating conductor in a rotating spacetime, it is useful to start with a simpler analogous 
configuration: that of a rotating magnetized conductor in a Minkowski (flat) spacetime. This will provide important insight for the search of 
general relativistic solutions and useful limits against which match the fully relativistic solutions. 

Consider therefore a conducting magnetized sphere of radius R rotating at angular velocity $7, and with the magnetic four-vector field 
B being uniform (in radius) inside the sphere and dipolar outside (This is a simple but instructive example.). Because of discontinuities in 
the fields across the surface of the sphere we will refer to as interior solutions those solutions valid within the radial range R IN < r < R, 
and to as exterior solutions those valid in the range R < r < oo. 



3.1.1 Interior Solution 

The interior solution for the electromagnetic fields of a magnetized sphere with magnetic moment aligned with the rotation axis was found 
by Ruffini and Treves in 1973 (Ruffini & Treves 1973). Extending it to the case of a misaligned rotator we obtain 

B T — (cos x cos + sin X sm $ cos A) , (48) 

B® — — (cos x sin & ~~ sm X cos 9 cos A) , (49) 
R d 

B^ = -|^sinxsinA, (50) 
R 6 

where /i is the magnetic dipole moment of the star, x ls the inclination angle of the magnetic moment relative to the rotation axis and 
\(t) = (f> ~ Qt is the instantaneous azimuthal position (see Fig. 1). 

The expressions for the components of the electric field are very simple to derive when one assumes that the sphere is a "perfect 
conductor" (i.e. a — > oo) and there are no conduction currents inside the sphere. In this case, Ohm's law can be used to obtain 

f firsinfl j 2^Slrsin6», . . . 

E = B = t (cos x sm " ~ sm X cos 6 cos A) , (51) 

c cR d 

_e f2rsin£>_f 2u£lrsm9 , . . „ . , 

hj = B = t (cos x cos o + sm \ sm " cos A) , (52) 

c cR d 

= . (53) 
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3.1.2 Exterior Solution 

The solution to this problem, i.e. to the form of the electromagnetic fields external to a misaligned rotating magnetized sphere, was found in 
1955 by Deutsch (Deutsch, 1955). The full solutions are complicated expressions involving spherical Bessel functions of the third kind, but 
they become much simpler when truncated at the lowest order. The magnetic field components, in particular, have the form 

B r = ~T ( cos ^ cos # + SU1 X sm & cos ^) ) (54) 
B® — (cos x sin 8 — sin x cos 6 cos A) , (55) 

= 4 sin x sin A , (56) 
while the corresponding electric fields are given by 

E r = — [cosx(3cos 2 6 — 1) + 3 sin x cos A sin 6> cos #] , (57) 



j2 



E = < — — cos x sin 8 cos 8 + sin x 

cr Z I r 2. 



?2 



■" " 1 t( c °s 8 — sin 9) 

r 2 



cos\} , (58) 



= ^sm X cos0^nx(l--%^ . (59) 

Three interesting features of solutions (54)-(56) and (57)-(59) should be noticed. The first one is given by the periodic time modulation 
introduced by the precession of the magnetic moment and which disappears when the dipole is aligned, i.e. for x = 0. The second feature 
is that, as one might have expected on the basis of symmetry considerations, the toroidal components of the external electromagnetic fields 
are just a by-product of the misalignment between the rotation axis and the magnetic dipole and again disappear when x = 0. Finally, the 
third relevant feature is the appearance of an electric field of O(0,) introduced by the rotation of the sphere and whose quadrupolar part [i.e. 
oc (3 cos 2 6 — 1)] is present also in the case of an aligned rotator. As we will see in Section 3.2.2, where we study the analogous problem in 
a slowly rotating spacetime, an additional contribution of 0(u>)^ to the form of the external electric field will be introduced by the general 
relativistic frame dragging effect § . 



3.2 Rotating Magnetized Conductor in a Slowly Rotating Spacetime 

We now consider the general relativistic analogue of the problem in Section 3.1 and look for a solution of Maxwell equations (30)-(33) and 
(34)-(37) assuming that magnetic field of the star is dipolar. To simply the search for a solution we look for separable solutions of Maxwell 
equations in the form 

B f (r,8,<p, X ,t) =F(r)*i(0,0,x,t) , (60) 
B § (r, 8, 0, x, t) = G(r)* 2 (e, 0, x , *) , (61) 

B*(r, 8, <p, x, t) = H(r)V 3 (e, <t>, X , t) , (62) 

where F(r), G(r), and H(r) will account for the relativistic corrections due to a curved background spacetime. 

A considerable simplification comes from the fact that, at first order in fi, the solutions for the electromagnetic fields will not acquire 
general relativistic corrections to their angular dependence. We therefore expect that, as for the case of the Deutsch solution, the general 
expressions for the angular eigenfunctions *i, with i = 1, . . . , 3, will have a complicated angular dependence expressed in terms of spherical 
Bessel functions of the third kind. This is however over-complicated and for most of the astrophysical applications it would sufficient to know 
the form at the lowest order which can be known by requiring that the solutions match the lowest order solution for a misaligned rotating 



t Hereafter we will refer to as 0(ui) any quantity that is the result of the dragging of reference frames and that is therefore oc <?04>- 

§ Note that an electric field induced by the rotation of the star must appear also in the general relativistic case. This is not present in the solution proposed by 
Prasanna and Gupta 1997, where the external electric field is only of 0(ui) and the radial dependence does not contain higher order terms in M/r. 
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dipole in flat spacetime. In this case, then, we obtain 

*i (9, (f>, x, t) = cos x cos 9 + sin x sin 9 cos A(t) 



(63) 



^2(9, <j>, x, t) = cos x sin 9 — sin x cos 9 cos A (t) 



(64) 



*3(0,X,*) = sinxsinA(t) , 
which also satisfy the following useful relations 

#1,9 = -*2 , #1,0 = -*3Shl0 , #2,9 = *1 , *2,0 = *3 COS ( 

Maxwell equations (30), (45)-(47) with the ansatz (60)-(62), yield the following set of equations 

^(r 2 -F) r + 2e A r(?J sm9 (cos x cos # + sin x sin 6< cos A) + e A r (H — G) sin x cos A = 0, 



(H — G) cos 6 sin x 
re*/J) + e* +A _F sin y sin x 



sin A 



cue 



47TCT 

,-* 



■ cos A 



= 4nr sin # J r 



we 

sin A cos A 

47TCT 



[£<»-°>l 



j—®,r (G — H) — ^ — 

47TO" 47TCT 



+ e* +A F 



cos 8 sin x sin A 



(re*(?) 



+ e* +A F 



(cos x sin S — sin x cos 6 cos A) 



-4ne rsm9J 



(65) 

(66) 

(67) 
(68) 
(69) 



(70) 



Next, we will distinguish between an external vacuum solution to Maxwell equations (for which fully analytic solution can be given) 
from the interior non-vacuum solution. Since we are treating the interior of the star as a perfect conductor and the exterior of the star as 
vacuum, we can impose J r = J = = in (67)-(70) and obtain as Maxwell equations for the radial part of the magnetic field 



(r 2 F) r + 2e A rG = 



+ e * +A F = . 



(71) 
(72) 



H- G = . 



(73) 



Note a first important result in the system of equations (71)-(73). In the case of stationary electromagnetic fields, the general relativistic 
frame dragging effect does not introduce a correction to the radial eigenfunctions of the magnetic fields. In other words, in the case of infinite 
conductivity and as far as the magnetic field is concerned, the study of Maxwell equations in a slow rotation metric provides no additional 
information with respect to a non-rotating metric. The frame dragging effects are therefore expected to appear at 0(lj 2 ). 



3.2.1 Interior solution 

Limiting the solution to an inner radius R IN removes the problem of suitable boundary conditions for r — > 0, and reflects the basic ignorance 
of the properties of magnetic fields in the interior regions of neutron stars. 

It is important to notice how the system of equations (71)— (73) combines information about the structure and physics of the star (through 
the metric functions $ and A) with information about the microphysics of the magnetic field (through the radial eigenfunctions F and G). 
As a result, a relativistic solution for the interior electromagnetic field cannot be given independently of a self-consistent solution of Einstein 
equations for the structure of the star. In practice, to calculate a generic solution to (71)-(73), it is necessary to start with a (realistic) equation 
of state and obtain a full solution of the relativistic star. Once the latter is known, the system of equations (71)— (73) can be solved for a 
magnetic field which is consistent with the star's structure and corresponds to a magnetic configuration of some astrophysical interest. (This 
is what done, for instance, by Gupta et al. 1998 in the case of an internal dipolar magnetic field). 

Alternatively, one might specify a magnetic field configuration and look for a compatible equation of state for the stellar structure (This 
is a less satisfactory way to proceed but one which is useful to get insight into the problem.). In this case, the simplest possible solution to 
the system (71)-(73) is one in which the magnetic field is constant throughout the region of the star of interest and is therefore the general 
relativistic analogue of the solution presented in 3.1.1. In this case, then 

F =§<", G = -^V = -e- A F, (74) 
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where Ci is an arbitrary constant whose value can be determined after imposing the continuity across the star surface of B r . 

We can now check whether the solution (74) is physically possible. Using (74) in the system of equations (71)— (73) requires that the 
metric functions satisfy the condition 

<I>+A 



[re ) 

Recalling now that Einstein equations for a spherical star yield 



0. 



e 2A W = l 1 _ 2m(r; 



r < R , 



with m(r) = 47r J q r 2 p(r)dr, and p(r) being the total energy density, we can rewrite (75) as 



2A / 

,r = e I 



m + rm r 



(75) 



(76) 



(77) 



On other hand, the solution of the Einstein equations for the interior of a relativistic spherical star (i.e. the solution of the Tolmann- 
Oppenheimer-Volkoff equations; Tolmann, 1939; Oppenheimer & Volkoff, 1939) requires that 



where P is the isotropic pressure. The comparison of (78) with (77) shows that the general relativistic uniform magnetic field solution 
1 



B r = 



(cos x cos 8 + sin \ sin cos A) C\p . 



B = - 



R 3 



B 4 



R 3 



(cos x sin # — sin x cos 8 cos A) C\p . 



(sin x sin A)Ci/i . 



(78) 

(79) 
(80) 
(81) 



is possible only for the (unrealistic) case of the "stiff matter" equation of state P = p. 

The corresponding form of the internal electric field is also straightforward to derive in the case of no conduction currents. In this case, 
in fact, Ohm's law (39) and (40) yield the simple expressions 



f _ ujr sm V g _ _ 



ci? 3 



-w (cos x sin 9 — sin x cos 8 cos A) Cip, , 



E § _ u>r sm V f 



e r sm ( 



cR 3 



■u) (cos x cos 8 + sin x sin 8 cos A) Cip , 



(82) 
(83) 



Ef = , (84) 

where we have taken into account that p e — O(co) and that the contribution proportional to Ljp e is therefore of higher order^. Note that, 
apart for red-shift correction proportional to e~*, equations (82)-(84) are the same as (51)-(53) with O, being replaced by the effective fluid 
velocity measured by a free falling observer ui. The internal charge density corresponding to the electrical field (82)-(84) can be calculated 
after imposing that [cf. eq. (34)] 

-A 



Using now expressions (82)-(84), we easily obtain 

-A 



3e 



-#- e ( 



, r rsinl 



sin0 I ■ a a w 2e - Kn 
^ 3 (cos x sm 8 — sm x cos 8 cos A) — w cos x }Cip . 



(85) 



(86) 



3.2.2 Exterior Solution 

The exterior solution for the magnetic field is simplified by the knowledge of explicit analytic expressions for the metric functions $ and 
A. In particular, after defining N = e* = e~ A = (1 - 2M/r) 1/2 , the system (71)-(73) can be written as a single, second-order ordinary 
differential equation for the unknown function F 

d_ 

dr 



2F = . 



(87) 



^ Prasanna and Gupta (1997) have used the assumption of infinite conductivity also for the matter outside the neutron star. We note that their expressions for 
the electric fields do not contain the (important) contribution of O(Q) and the radial component does not seem to satisfy Ohm's law. 
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Introducing now the new variable x = 1 — r /M, equation (87) can be written as 



dx 



(88) 



1 — x J dx 

Equation (88) is an example from a class of equations which can be solved in terms of the Legendre functions of the second kind Qi (see 
Appendix B for details on the derivation of the solution). In the case of equation (88) we have t = 1, and (Jeffrey 1995) 



Qi = g ln 

The radial eigenfunctions F(r), G(r), and H(r), are then given by 

3 



1 +x 
1-x 



1 . 



F(r) = - 
G(r) = 



4M 3 
3A 



, ,r2 2M 

ln A H 



1 + ^ 

r 



4M 2 r 



i^ lniV2 + ^ + 1 



ff(r) = G(r) , 

and satisfy the following boundary conditions: 

(i) refer to a vanishing field at infinity, i.e. 

lim F(r) = , 

r — ► oo 

(ii) reduce to a flat spacetime solution for a dipole, i.e. 

lim F{r) = %, 

M/r^0 r 3 



lim G(r) = ; 



lim G(r) 

M/r->0 



(89) 

(90) 

(91) 
(92) 

(93) 
(94) 



(Hi) coincide with the corresponding radial eigenfunctions found for a Schwarzschild spacetime (Ginzburg & Ozernoy, 1964; Anderson 
& Cohen, 1970). This is what we expected since there are no first order contributions due to the rotation of the spacetime. 
Using expressions (90)-(92) we can now determine the value of the matching constant Ci by requiring that the radial magnetic field is 
continuous across the star surface, i.e. that [B r (r = R)] IN = [B r (r = R)] EXT . As a result, we obtain 



Ci = - 



3iT 

4M 3 



In 



, 2M\ 2M ( 
1-^ +^1 



R J 



R 



+ 



M 



F(R)R 3 



whose flat spacetime limit is 



lim Ci = 2 . 

m/r->o 



(95) 



(96) 



Collecting all the expressions for the radial eigenfunctions, the stationary vacuum magnetic field external to a misaligned magnetized 
relativistic star is given by 



B r 



4M 3 

3N 
4M 2 r 

3N 
4M 2 r 



2M 



ln N + 1 + — 



M 



(cos x cos 6 + sin \ sin cos A) fi , 



— IniV 2 + -t- + 1 



M 



N 2 



(cos x sin 9 ~ sin x cos 9 cos A) n . 



(sin x sin A)/i . 



(97) 
(98) 
(99) 



The search for the form of the electric field is much more involved than for the magnetic field. However, hereafter we will make use of 
the insight gained in Section 3.1.2 as a guide and start the derivation of the solution by rewriting vacuum Maxwell equations (31)— (33) and 
(34) as 



3u>r 
4M 3 iV 
3uir 



4M 2 



, * r 2 2M 
ln N 2 H 



1 + ^ 
r 



9ujr 
4M 3 ' 



, »,2 2M 
InN 2 H 



1 + ^ 

r 



3lj 



4M 2 



sin x sin 6 sin A 
sin x sin # cos S sin A 
(cos x cos 9 + sin x sin cos A) sin # 
sin x cos A 



(sintf^) 



-E" 



r i 1 
— IniV 2 + -2+1 
M iV 2 



A 



sin (r 2 £") + r (sin 9E § ^j ^ + rE^ 



E r 



rNE 
0, 



(100) 
(101) 

(102) 
(103) 



and which already indicate that the dragging of inertial frames with angular velocity lo introduces electric fields in the surrounding space 
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when magnetic fields are present. Using as a reference the solutions (57), (58), and (59) for a misaligned rotating sphere in Minkowski 
spacetime, we look for the simplest solutions of vacuum Maxwell equations in the form 



E r = (/i + /s) cosx(3cos — 1) + (gi + gs) 3 sin \ cos A sin 8 cos 8 , 



(104) 



E 6 — (/ 2 + fi) cos x sin 9 cos 8 + (c? 2 + gi) sin \ cos A — (g 5 + go) (cos 2 8 — sin 2 8) sin \ cos A , 



(105) 



E = [95 +96- (92 + gi)] sin \ cos 8 sin A , 



(106) 



where the unknown eigenfunctions /i — /4, and gi — go can be found as solutions to vacuum Maxwell equations and have radial dependence 
only. Substituting (104)-(106) in (100)~(102) we obtain the following set of linear differential equations 



N (r 2 /i) r + r/a = ; 



(107) 



(rJV/ a ), r +6/i = ; 



(108) 



N (r 2 / 3 ) +rf 4 =0. 



(rJV/ 4 ), r + 6/ 3 -^ 



, Ar2 2M / M 
In N 2 H 1H 



= . 



N (r 2 5l ) r + 2rg 5 = ; 



(109) 
(110) 
(111) 



{rNg 5 ) r + 3 gi = 0, 



(112) 



N (r 2 53 ) r + 2rg 6 = ; 



(113) 



(riV fl6 ) r + 3 5 3-^3M 



, , r2 2M / M 
In 7Y H 1H 



= . 



(114) 



Note that both the sets of radial eigenfunctions fi — f^, and gi, gz, gs, ge are linearly independent, but that relations can be written between 
the two sets. In particular, the comparison of equation (107) with (111) and of equation (109) with (113) indicates that 



9i= h, 



93 = h 



55 = 



h 



56 = 



(115) 



We start the search for explicit expressions for the radial eigenfunctions by combining equations (107) and (108) to obtain a single 
differential equation of second order for the unknown function /i 



dr 



[('-^)|C«]-« 



. 



(116) 



and which can again be recast in a form similar to equation (88). Proceeding in a way analogous to what done for the magnetic field (see 
Appendix B for details) it is possible to realize that the solution should be expressed in terms of a Legendre function of the second kind and 
of order 1 = 2. Recalling now that (Jeffrey 1995) 



Q 2 (,) = I(3, 2 -l)ln(j±j 



3a; 
T ' 



we obtain, as solution to (116) at the i = 2 order in the expansion 



h 



6cR 2 



C1C2 



„ 2r\, Ar2 2M* 2M A 
Ml 3r 2 r 



(117) 



(118) 



where C 2 is an arbitrary constant to be determined through the imposition of boundary conditions. Making now use of the equation (107) we 
also obtain that 



h = -^C,C>N 



IniV 2 



2M/_ 

3r 2 N 2 



fi . 



(119) 
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In a similar way, the solutions to equations (109) and (110) are found to be 



h = 



u = - 



8M 5 c 



Ca 



NIC, 



2r 
M 



In N 



2M 2 
3r 2 



In AT -2- 



2M 
r 


-4 


2M 2 
5r 2 


2M 2 ' 


+ 


4M 4 


3r 2 N 2 


15r 4 iV 2 



lniV 2 + — — 

5r J 



(120) 
(121) 



where again C3 is an arbitrary constant determined through the boundary conditions. Finally, the functions gz and qa are given by 



92 



3fir 

8M 3 cN 



, 2M 
lniV 2 H 



1 + ^ 
r 



ff4 = ~n 52 = ~8M^ 



In iV 2 H 



1 + 



M 



(122) 
(123) 



so that 32 + 54 = (w/fi)<?2- 

Collecting again all the expressions for the radial eigenfunctions, the stationary vacuum electric field external to a misaligned magnetized 
relativistic star is given by 



E r = 



15^ 
16M 5 c 

n 



c 3 



3- 



2r 



, , r2 2M 2 2M A 
In Af 2 + — - H 4 



C1C2 



6ci? 2 
I 45wr 3 



3 " S) 1,1 V "- 

M 



3r 2 



2 2Af 2 2M ; 
- ■ 1 4 

3r 2 r 



2M 2 , „ r2 4M d 



"| 16M 5 c 



+ 



2ci? 2 

3u>r 
8M 3 cN 



C1C2N 



\ Ml 
V M) 



In Af 2 - 2 



In A> -2- 



3r 2 Af 2 
2M 2 



In N 2 H 



1 + ^ 
r 



3r 2 AT 2 
(sin x cos A)/i 



[cos x(3 cos 2 5 — 1) + 3 sin x cos A sin 6 cos 0] fi , 
4M 4 1 

[2 cos x sin 9 cos — (cos 2 6 — sin 2 0) sin x cos A] 



(124) 



(125) 



45wr J 
16M 5 c 

3wr 
~8M 3 cN 



N{C 3 



lnA/ 2 -2- 2M 



lniV 2 + ^|l + 



M 



3r 2 A T2 

(sin x cos # sin A)p 



4M " 4 1 ^ ^ ^ »r 



lnA/ 2 -2- 2M 



3r 2 iV 2 



(126) 



As anticipated in Section 3.1.2, expressions (124)-(126) confirm that the general relativistic dragging of reference frames introduces a 
new contribution to the form of the electric field which does not have a flat spacetime analogue. This effect is 0(uj) and therefore present 
already in a slow rotation approximation. This is in contrast with what happens for the magnetic fields, where higher order approximations 
of the form of the metric are necessary for frame dragging corrections to appear. 

The values of the arbitrary constants C2 and C3 can now be found after imposing the continuity of the tangential electric field across 
the star surface. Using then (82)-(84) as solutions for the internal electric field and imposing that [E 6 (r = R)] IN = [E 6 (r = R)] EXT as 
well as [E*{r = R)] IN = [E*(r = R)] EXT , yields 



C 2 



1 

Af 



R_ 

M 



lnAT^-2- 



_2W_ 

3R 2 N 2 



C 3 



2M 2 
15i? 2 



, » r 2 2M 



(127) 
(128) 



with N 2 = N 2 (r = R) = 1 — 2M/R. It is now also possible to calculate the surface charge distribution <j s resulting from the discontinuity 
across the star's surface of the radial electric field. Explicit expressions for this, as well as for the surface currents corresponding to the 
discontinuities across the surface of B e and B^, will not be given here but can be found in Appendix C. 

Before concluding this Section on stationary solutions we will comment on the relevant limits of equations (124)-(126). Firstly, we 
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verify that they reduce to the Deutsch solutions (57)-(59) in the limit uo — and M/r, M/R —> 0. In this case, in fact 



lim /i (r) 

M/r, M/R—>0 



lim f 2 (r) 

M/r, M/R->Q 



lim g 2 (r) 

M/r, M/R^O 



lim g 4 (r) 

M/r, M/R—>0 



M/r, M/R^O 



= lim Qi(r) = lim 

Cr* M/r, M/K->0 M/r, M/R->0 



2/ittR 2 



lim / 3 (r) = 



lim 

M/r, M/fl^O 



fi(r) = lim g 3 (r) = lim g 6 (r) 

M/r, M/R-^Q M/r, M/R-,0 



Secondly, in the limit of an aligned dipole in a Schwarzschild spacetime, x = = w, and equations (124)-(126) reduce to 



(129) 
(130) 
(131) 
(132) 
(133) 



Q.R 



, ,r2 2M 

l « N » + -R 



1 + 



M 



1 - 



R_ 

M 

2r 



3^i? 
4M 3 ciV2 



, »r2 2M 



1 + 



3 -M llniV + 



2 2M Z 2M ' 
- ' — _ -| 4 

3r 2 r 



(3 cos 2 6 — , 



1-^-7 IniV 



Af 



2M 



2 -i -i 



V Ml 



In AT 



3R 2 Nl, 



2M 2 



3r 2 N 2 



(sin^cos^)^ , 



(134) 



(135) 



: . 



(136) 



Note that (134)— (136) do not coincide with the corresponding expressions found by Sengupta (1995). A straightforward calculation would 
show that his suggested expressions, while reducing to the Deutsch solution in the flat spacetime limit, do not satisfy Maxwell equations. A 
possible explanation for the disagreement could be found in the method followed by Sengupta in his derivation which is not based on the 
explicit solution of Maxwell equations. Because of this, subsequent results obtained on the basis of Sengupta's expressions for the external 
electric field (e.g. De Paolis et al., 1999) should be revisited in terms of expressions (134)— (136). 



4 NON-STATIONARY SOLUTIONS TO MAXWELL EQUATIONS 



In this Section we will drop the assumption of infinite conductivity which prevented the variation of the star's magnetic moment and led to 
the stationary electromagnetic fields presented in the previous Sections. Here, on the contrary, we are interested in time evolving electromag- 
netic fields and, in particular, in establishing the general relativistic corrections to the induction equation. A direct consequence of a finite 
conductivity is, in fact, the generation of a time varying charge density and conduction currents which will be then responsible for the Ohmic 
decay. Using Maxwell equations (34) and Ohm's laws (39)-(42), we find that the space charge density p e = p e {t, r, 6, <f>) inside the star has 
a zeroth-order contribution given by 



Pe = 



16n 2 ar 2 sin6> 



(sinOB^ 



rB 



+ 



,6 

e B j. — sin t 



-* i 
e + 



(e*r£«) 



*+A R f 

e Bo 



ce A / cos 9 sin x sin A \ (G — H 



16tt 2 c 



< } +0(SI) ; (137) 
(138) 



where the second expression is the one obtained after using the ansatz (60)-(62). It follows from equations (143) and (144) that the zero-order 
term in equation (138) vanishes, so that the leading contribution is at first order in ft. Using now equations (31), (46), (47), (137) and Ohm's 
laws (41), (42), we obtain the evolution equation for the radial component of magnetic field 
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dB r 



1 



Anar 2 sin 6 ) sin 8 



S+A or 



(e*rB*) -e* +A iT 6 







ujr 






sin 6* 




[( 






Ana 





47TCT 



*+A 7-,f 

e i) a sin 



+ < sin t 



(e*rB*) 
(e*rB*) 



47Tcr sin 9 

sinOB^ -B § 
I .e 

fie"* ' 



B^-sin^fe*^ 



Ana 



■(smOB^ 



rB 







sin<9 


.... 







(e*rB 5 ) 



(139) 



Similarly, using equations (32), (45), (47), (137) and Ohm's laws (40), (42), we obtain the evolution equation for the polar component 
of magnetic field 



dB e 
dt 



Anar 



(e*rB § ) 



$+A R f 

e Be 



+ - <; uje 

Ana 



-(*+A) 



e* +A B r 4 — sin^ 



]}„ + { 



Ana 



B 



Ana 



(sin 0-5"*) 



rB 



+ e 



-(*+A) 



e #+A B f 



sinfl (eVe* 



+ <e 



Anar 2 sin 2 6 1 



(sin 6^) 



B 



-(*+A) 



47T(T 



(sin 6^) 



*+A R f 



QB" 



(140) 



Finally, equations (33), (45), (46) and the Ohm's law (40), (41) yield the evolution equation for the toroidal part of the magnetic field 



OB* 



Anar sin 6 



#+A r>f 





f we -(*+ A ) r 


L-l 


47TCT 







(1 


sine 



(smOB^ 



B 



Ana sin 



(sin^) 



*+A D f 



73" 



< (re*B*) 



Oe" A sin 6> 



(r 2 /3 f ) +fi (sin^) 



(141) 



A first important feature of equations (139)— (141) is that besides the relativistic corrections due the monopolar part of the gravitational 
field (proportional to M/R and already present in the non-rotating case), the rotation of spacetime introduces additional corrections related 
to the dipolar part of the gravitational field (and proportional to cu) to the decay of the magnetic field. A second relevant aspect of equations 
(139)— (141) is that they do not show the degeneracy encountered in the time evolution of the magnetic field in a non-rotating spacetime [cf. 
equation (146)]. In that case, in fact, the three induction equations for the components of the magnetic field reduce to a single evolution 
equation (Sengupta, 1997; Geppert et al. 2000). Here, instead, the three equations remain distinct and a particular field component might be 
favoured during the decay. Finally, equations (139)— (141) do not factor out the angular part as it is the case for a non-rotating, aligned dipole 
[cf. equation (146)] and the evolution of the magnetic field has therefore properties which depend on the angular position in the star. As a 
consequence of this, an initially dipolar magnetic field might not remain as such during its decay. This could be relevant for the evolution of 
the magnetic field of pulsars and more particularly of magnetars. 
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Using now the ansatz (60)-(62), we can write equations (139)— (141) in the more compact form 

^e*r (G — if)J sin x cos A 



dF c 2 e' A 

(cos 8 cos v + sin ft sin y cos A) sin 9 = 

Ot ' Altar 2 



( C *rG) 



+ e*+ A F 



sin 6 (cos 9 cos x + sin ft sin x cos A) 



--!-sinxsinA{[a;r( J ff-G)] r (l-2sin 2 6»)+2we-* (e*r#) + e* +A F 
47rcr |^ [V /,r 

-fir (G $ r (l-2sin 2 ft) jj , 



(142) 



dG 
dt 



(sin # cos x — cos # sin x cos A) 



e* (G - H) 



cos ft sin x 



(e*rG) 



+ e*F 



47rar ^ r sin 2 ft 

(sin ft cos x — cos 9 sin x cos A) 



cos A H — - sin A 



47TCT 



47ra 



„2„-A 

■ sin A = 

dt i-Kcrr 



cos ft sin x sin A< e [wr (G — H)] r + w 



F + e -( A +*) (re*tf) + fi [$ r e- A r (G - tf) 



(«**■*), 



+ e*F 



we 

sin A cos A 

47T(7 



+ 



c 2 e f (G - H) 
4-Kar 2 sin 2 9 



uie 

sin A cos A 

47TCT 



(143) 
(144) 



where F, G, and H satisfy the constraint equation (30) 

[(r 2 F) r + 2e A rG] sin ft (cos x cos ft + sin \ sin ^ cos A) + e A r (H - G) sin x cos A = , (145) 

The set of equations (142)-(144) is too complicated to be solved analytically even when analytic expressions are available for the metric 
functions (e.g. for a constant density stellar model). The numerical solution of (142)-(144) for a number of equations of state together 
with a self-consistent evolution of the star's angular velocity and electrical conductivity will be presented in a separate paper (Rezzolla et 
al. 2000). Note that equations (142)-(144) could also be derived through a vector potential defined so that the electromagnetic tensor 
F M „ = A v ^ — A^.y. Details of this derivation can be found in Appendix D. 

An interesting limit of the induction equations (142), (143) is the one for a non-rotating dipole in a spherically symmetric spacetime. In 
this case: Q. = = u>, x — 0, and H, H/dt are not determined [cf. equations (62) and (65)]. As mentioned before, the induction equations 
are degenerate in this case and the unique evolution equation is then 



dF = c 2 e~ A 
dt 4nar 2 



2e* +A F 



(146) 



corresponding to the solution found by Geppert et al. (2000). When the metric functions $ and A refer only to the vacuum region of spacetime 
external to the star, equation (146) further simplifies to 



OF _ c 2 
dt Anar 2 



2M 



2F 



(147) 



and which now corresponds to the solution found by Sengupta (1997). 



5 CONCLUSION 

We have presented analytic general relativistic expressions for the electromagnetic fields internal and external to a slowly-rotating magnetized 
neutron star. The star is considered isolated and in vacuum, but no special assumption is made on the orientation of the dipolar magnetic 
field with respect to the rotation axis. The solutions to Maxwell equations have been considered both for an infinite and for a finite electrical 
conductivity. 

In the first case, corresponding to stationary magnetic fields, we have shown that the general relativistic corrections due to the dragging 
of reference frames are not present in the form of the magnetic fields but emerge only in the form of the electric fields. In particular, we 
have shown that the frame-dragging provides an additional induced electric field which is analogous to the one introduced by the rotation of 
the star in the flat spacetime limit. In the case of finite electrical conductivity, on the other hand, corresponding to decaying magnetic fields, 
we have shown that corrections due both to the spacetime curvature and to the dragging of reference frames can be found in the induction 
equation. An interesting result obtained in this regime is that the rotation of the star eliminates the degeneracy in the components of the 
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induction equation which remain therefore distinct. Furthermore, rotation and dipole misalignment do not eliminate the angular dependence 
in the induction equation and, as a result, an initially dipolar magnetic field might evolve towards a different configuration during its decay. 

Because of their complexity, the evolution equations found for the magnetic field require a numerical integration which will discuss in 
detail in a forthcoming work (Rezzolla et al. 2000). There, we will also present direct comparisons between the flat and the curved spacetime 
solutions and quantify more precisely the importance of the general relativistic corrections. 

One of the relevant aspects of the solutions presented in this paper is that they provide a lowest order analytic form for the electromag- 
netic field in the spacetime of a slowly rotating misaligned dipole subject to assumptions which, while giving simplifications, allow the major 
features of a realistic solution to be seen. In this sense, they reflect a rather general physical configuration and could therefore be used in a 
variety of astrophysical situations. 
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APPENDIX A: THE ELECTROMAGNETIC TENSOR 

For completeness, we provide below the explicit expressions for the components of the electromagnetic tensor used throughout the paper. In 
a coordinate basis, and at first order in Q., these components are given by 

/ -e*E r -iue A r 2 sin 9B e -e*E e + weV sin 6B r -e* ' E^ \ 



"a/3 



e*B r +weV smOB 6 
e* E e - Lue A r 2 sm6B r 







e A r 2 sin6»B e 



-e A r 2 sm6B r 



-eV sines 6 
e A r 2 sin 0B r 




V 



(Al) 



/ 



The matrix (Al) can also be expressed in terms of the electromagnetic field measured by the ZAMO observers, in which case it takes 
the form 



F af) = 



e 9+A E f +u)e A rsmeB e 
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Finally, we note that the components of the electromagnetic tensor in the ZAMO frame can be derived from (Al) with the transformation 



to obtain 



F - — pV-p" F 

/ -cE f -cE § 
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APPENDIX B: RADIAL EIGENFUNCTIONS 

In this Appendix we briefly sketch the procedure for the calculation of the radial eigenfunctions and that have lead to the solutions (90)-(92), 
(1 18)— (1 19), (122)-(123). In general, we look for a solution of the equation 



s{(i?f)sK 1 -*>'"]} + '" +1 »— 



for the function qe. Equation (Bl) can also be written as 

(1 - x 2 ) 4 - 2(1 + 2x)qt + [£ (£ + 1) - 2] q e = , 
with the dash representing a total derivative with respect to x. The solution of (Bl) has then form 



qe = 



d_ 

dx 



where Qe are Legendre functions of second kind (Jeffrey 1995). A proof of this comes from substituting (B3) into (B2) to obtain 
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It is now easy to realize that (B5) is identically satisfied since the content of the square brackets is, in fact, the second derivative of 
Legendre's equation 

(1-X 2 )q' e ' -2xq' e + £(£+ l)q t = 0, (B6) 
while all the remainder of (B5) is the first derivative of Legendre's equation (B6). 



APPENDIX C: SURFACE CHARGES AND CURRENTS 

We here give explicit expressions for the surface charge distribution a s and surface currents i k resulting from the discontinuities across the 
star's surface of the r— component of the electric field and of 6, 0— components of the magnetic field. 

Defining now [A]^ = [A(r = R)] EXT — [A(r = R)] IN , the surface charge density a s can be found as 
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where a;/; = uj(r = i?) and lj r = oj(r = R) 
In a similar way, imposing that 
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APPENDIX D: AN ALTERNATIVE DERIVATION OF THE INDUCTION EQUATION 



To confirm the results presented in Section 4 and to compare with the results presented in the literature in the case of Schwarzschild back- 
ground spacetime (Sengupta, 1997) we here present a derivation of equations (142), (143), (144) in terms of a vector potential A a defined 

as 



F a a = Ap, a - A a . 



(Dl) 



The use of a vector potential is sometimes looked at with skepticism (Geppert, Page and Zannias, 2000) in view of the non-commutativity of 
the covariant derivative, which could lead to ambiguities if Maxwell equations are expressed through double covariant derivatives of a vector 
potential. All of these ambiguities, however, are easily removed if the vector potential is introduced only when Maxwell equations (15) are 
recast in a form not involving covariant derivatives 
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In our derivation we start by using equations (7), (15), (17), neglecting the displacement current and taking the four-velocity of the 



conductor 11° in the form (38) [i.e. neglecting terms proportional to g F ob u w 0{u> )] we get 
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Using now (D2) we can write the general expression for the evolution of the vector potential Ai as 
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Using (D5), the induction equation for the evolution of the (f> component of vector potential can be written as 
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which coincides with equation (11) of Sengupta's 1998 paper when \ = = ui and e 2 * = N 2 = 1 — 2M/r. Similarly, the evolution 
equations for the other components of the vector potential are given by 
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In the case of a misaligned rotator, the explicit expressions for the "magnetic" components of the electromagnetic tensor are 

F r e = e A rH sin x sin A , 
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F, $4, — Fr sin 6 (cos 6 cos x + sin \ cos A sin 8) , 
using which the induction equations (D6)-(D8) assume the form 
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Using now the definition (Dl) it is possible to show that equations (D12), (D13), (D14) are equivalent to equations (142), (143) and 
(144) derived in the main text. 
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